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Abstract. 

We study the transport properties of a long non-uniform quantum wire where the 
electron-electron interactions and the density vary smoothly at large length scales. We 
show that these inhomogeneities lead to a finite resistivity of the wire, due to a weak 
violation of momentum conservation in the collisions between electrons. Estimating 
the rate of change of momentum associated with non-momentum-conserving scattering 
processes, we derive the expression for the resistivity of the wire in the regime of weakly 
interacting electrons and find a contribution linear in temperature for a broad range of 
temperatures below the Fermi energy. By estimating the energy dissipated throughout 
the wire by low-energy excitations, we then develop a different method for deriving the 
resistivity of the wire, which can be combined with the bosonization formalism. This 
allows us to compare our results with previous works relying on an extension of the 
Tomonaga-Luttinger model to inhomogeneous systems. 
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1. Introduction 

Recent experiments on quantum wires and carbon nanotubes PjElEllHElEllTllHlEl 
[TOl [TTl [T2| IT3| [HI |15] have stimulated a lot of interest in the transport properties of 
one-dimensional conductors. From a theoretical point of view, interacting electrons in 
one dimension form the so-called Luttinger liquid [16l[I7j, whose properties qualitatively 
differ from the conventional Fermi liquid state. Recent progress in fabrication techniques 
has made possible the experimental observation of various characteristic signatures 
of the Luttinger liquid, such as the power-law behavior of the tunneling density of 
states [121 lEl El, or the existence of separate spin and charge excitations p^. It 
is also expected, within the Luttinger liquid theory, that the dc conductance of a 
quantum wire connected to Fermi liquid leads is given by the quantum of conductance 
Go = 2e^//i [m [ini Eni- This quantization of the conductance has been reported in 
various experimental setups since its first observation in a quantum point contact [Il[2]. 

However, in a number of recent experiments [31 HI El El [TJ El [HI [TOl [H] , significant 
deviations from perfect quantization have been observed in the regime of low electron 
density. These deviations take the form of a shoulder-like structure below the first 
plateau of conductance. Although weak at the lowest temperatures available, this 
feature becomes more significant as the temperature increases, turning into a quasi- 
plateau at about 0.7 x {2e^/h). This so-called "0.7 structure", which is not expected 
in the Luttinger liquid theory, generated much theoretical interest, though there is at 
present no generally accepted microscopic theory. Most commonly, the experimental 
results are interpreted as originating from a spin-dependent mechanism. Such scenarios 
rely on a spontaneous spin polarization of the wire [3l [211 [22] , or on the existence of a 
local spin-degenerate quasi-bound state whose screening would lead to Kondo-like effects 
[23| [2^. Other proposals considered various scattering mechanisms involving plasmons 
[25] , spin waves [26] or phonons p7| . Several authors have also suggested that electron- 
electron interactions may affect the transport properties in quantum wire devices in a 
way that would be consistent with the "0.7 structure" [28 l l29l [30l [3T]. 

In this context, a number of recent theory papers studied the electronic transport 
in a quantum wire modeled as a one-dimensional system in which the interactions are 
limited to a small region between two non-interacting leads. They concluded that 
the backscattering of either single electrons or pairs were the only mechanisms to 
significantly affect the transport properties of the system [301 EI]; but only if the size 
of the interacting region is comparable to the Fermi wavelength of the electrons in the 
wire. If, on the other hand, the interaction strength varies smoothly over a much larger 
distance, such backscattering processes only lead to exponentially small contributions 
which can be neglected. Using the model of a non-uniform Luttinger liquid with position- 
dependent parameters, it was found that no correction to the quantized conductance of 
the wire arises in this regime [TH [TOl [20] . 

In this paper we show that even when the backscattering processes can be ignored, 
the non-uniformity of the interaction potential throughout the wire leads to a finite 
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resistivity at non-zero temperatures. Indeed, the inhomogeneity of the interaction 
potential breaks the translational invariance of the system, allowing for two-particle 
scattering processes that conserve energy but not momentum. In section |2l we 
qualitatively show how some of these processes give rise to a finite resistivity and perform 
the corresponding calculation in section [3l In section HI we present an alternative 
derivation of the resistivity in the language of the inhomogeneous Luttinger liquid 
model, allowing us to compare our results with previous works relying on this formalism 
[T8l flQl [20] . Finally, in section [5] we discuss the relation of our results to the experiments 
probing the transport properties of inhomogeneous quantum wires. A brief summary of 
some of our results was reported in [32] . 

2. Qualitative picture 

Let us consider an infinite one-dimensional system of weakly interacting electrons with 
a quadratic dispersion = /2m. To develop a qualitative picture of the physics 
involved, we restrict ourselves to the simple model of spinless electrons, with a uniform 
density n throughout the device. (We will tackle more realistic systems in the next 
section.) The inhomogeneity of the system comes from the electron-electron interaction 
whose strength varies smoothly along the wire. 

When one enforces a dc current / to flow through the device, the electrons start 
moving and acquire a drift velocity Vd proportional to this applied current: Vd = I/ne. 
In the reference frame moving with velocity Vd along the wire, the electronic subsystem 
is in an equilibrium state characterized by a Fermi energy ep and a temperature T. 
This was recently pointed out [33] in the context of Coulomb drag between two parallel 
wires. 

As we are interested in the low-energy properties of the system, we focus on 
temperatures T -C e^?, so that the only relevant excitations are close to the Fermi level. 
As a result, one can isolate two well-defined branches corresponding to two species of 
fermions: the right- and left-moving electrons. Within each branch, the velocity of the 
electrons can be approximated by a constant and is given by +vf and —vp respectively 
for right- and left-movers. Upon changing from the moving to the stationary frame of 
reference, the electron velocities are modified in order to account for the drift velocity, 
and change from ±vf to ±VF+Vd. The consequences for the electron fluid as described in 
the stationary frame of reference are two-fold. First, we need to introduce different Fermi 
energies for right- and left-moving electrons, ep — * e^'^ = (l/2)m(fi? ± Vd)"^. Second, 
since the density of states at the Fermi level is inversely proportional to velocity, we now 
have different densities of states for the two subsystems, u oc l/vp ^ ur^l oc 1/ {vF^Vd)- 

The latter result implies that the energy spacing between states is not only modified 
as we change the frame of reference, but also differs between the right and left branches 
in the stationary frame. Compared to the moving frame, the energy levels are stretched 
near the right Fermi point. This results in a somewhat broader distribution function, 
which can be interpreted as a slightly higher effective temperature Tr for the right- 
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Figure 1. Electronic spectrum in tlie stationary frame, with the Fermi energies 
and effective temperatures for the right- and left-moving branches. The corresponding 
distribution functions near the right and left Fermi points are displayed as functions 
of energy. An example of non-momentum-conserving scattering process is provided. 



moving electrons (see figure [2]). Similarly, near the left Fermi point, the energy levels are 
squeezed compared to the moving frame, resulting in a narrower distribution function, 
corresponding to a lower effective temperature Ti for the left-moving electrons. These 
effective temperatures follow the change in the density of states and are given by: 



The nature of these effective temperatures can be understood formally, by noticing that 
in the stationary frame, the system is no longer in thermal equilibrium because of the 
finite electric current. It follows that, quite generally, the occupation probability of a 
given state is no longer given by the standard Fermi-Dirac distribution. However, the 
introduction of the effective temperatures ([1]) for right- and left-movers enables one to 
write their occupation probabilities as Fermi functions of energy. 

Because right- and left-movers have different temperatures, it is natural to 
expect that electron-electron interactions will give rise to thermalization between the 
two branches. In a uniform system, two-particle scattering processes cannot lead 
to thermalization as the conservation of both energy and momentum only allows 
processes which either exchange the momenta of the two incoming electrons or leave 
them unchanged [M]. On the other hand, in the case of inhomogeneous wires, the 
strength of the interaction potential is non-uniform so that the system is no longer 
translationally invariant, and two-particle scattering processes which conserve energy 
but not momentum are allowed. 

A typical example of such electron-electron scattering processes is shown in figure [H 
It describes the scattering of two electrons from an initial state with momenta p and k 
to a final state with momenta p' and k', and violates the momentum conservation: 
p' + k' — p — k = P < 0. Though the loss of momentum associated with this 
scattering process may affect the transport properties of the system, one could argue 
that it is compensated by an equal gain of momentum corresponding to the inverse 




(1) 



Electronic transport in inhomogeneous quantum wires 



5 



process (p', k') — > (p, k). This is however not the case here because of the temperature 
difference between the two branches: the processes involving a transfer of energy from 
the "warmer" right-moving branch to the "colder" left-moving one statistically occur 
more often than the corresponding inverse processes. As a result, the electronic system 
loses more momentum than it gains. 

This overall loss of momentum can be viewed as resulting from a damping force, 
associated with the electron-electron collisions, and proportional to the temperature 
difference between the right- and left-moving branches. In order for a constant current 
to flow through the wire, this damping force has to be compensated by a driving force. 
The latter originates from a local electric field which appears as a response of the 
system to the external current. Using the force balance, and keeping in mind that the 
temperature difference Tr — Tl oc /, this local electric field is proportional to the applied 
current bias. This implies a finite resistivity of the wire. 

3. Weakly interacting electrons in the stationary frame 

The above arguments provide a physical picture of how inhomogeneities lead to a finite 
resistivity. We now proceed with the calculation of the resistivity. 

3.1. Model 

Our starting point is a one-dimensional system of weakly interacting electrons with 
spins. In order to account for a non- uniform electron density n{x). we introduce a one- 
particle potential U{x) originating from the surrounding gates and impurities in the 
substrate. Moreover, the interaction between electrons is inhomogeneous, and described 
by a smoothly varying potential y(r, i?), given in the center-of-mass coordinates. The 
Hamiltonian for this system takes the form 

Ho + H,^t (2a) 

^0= 'Zjdxi^\{x)[-^ + U{x)-,^i^,{x) {2b) 

Hint = \Y.j^Rj^^ nr,R) ^\ (it:+ 4 (i?- 

where il^\{x) creates an electron with spin projection 7 at position x, and 11 is the 
chemical potential. We assume that the potential U {x) is a smooth function of position, 
and that U{x) <S This allows us to introduce a position-dependent Fermi energy 
ep{x) = — U{x). Similarly, the position-dependent Fermi momentum and velocity are 
straightforwardly defined as Pf{x) = \j2meF{x) and vf{x) = pF{x)/m. 

We keep a very general form for the interaction potential between electrons and 
only make the following assumptions concerning its characteristic length scales. On the 
one hand, we assume for simplicity that the interaction is short-range, the potential 
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Figure 2. The three non-momentum-conservmg processes that contribute to the 
resistivity. These processes can be designated using the standard notations involving 
coupling constants [35]: the scattering process represented in (a) corresponds to 
92\\ ~ while (b) corresponds to g2±, and (c) to gi±. Following this correspondence, 
in the text we use the notations ||, 2_L and 1_L to refer respectively to (a), (6) and (c). 



decaying rapidly as a function of the distance r between electrons. On the other hand, 
since we consider a non-uniform system, the interaction depends on the position R 
of the center of mass. The variations with respect to R are smooth and occur at a 
typical length scale d, large compared to both the Fermi wavelength and the range of 
the interaction potential. Similarly, we assume that the potential U{x) varies at the 
same typical length scale d as the interaction strength. 

3.2. Resistivity 

We focus now on temperatures in a broad range hvp/d ^ T ^ ep- In order to compute 
the resistivity of the wire, we consider a force balance on a small isolated segment of 
wire taken at position x, whose length Ax well exceeds the range of the interaction while 
satisfying hvp/T -C Ax <C d. When an external current / is applied to the device, the 
response of the system manifests itself as a local electric field E{x) = p{x)I, which in 
turn leads to a driving force eE{x)n{x)Ax acting on the electrons. This driving force is 
compensated by a damping force AF resulting from the inhomogeneous electron-electron 
interaction, so that the resistivity can be written as 



The damping force can be evaluated as the change in momentum per unit time 
associated with two-particle scattering processes. In the regime kpd ^ 1, the processes 
with a large momentum difference compared with the Fermi momentum lead to 
exponentially small contributions. As a result, in what follows we focus on processes 
which only weakly violate the momentum conservation (see figure [2]). 

Because of the non-uniformity of the wire, strictly speaking the momentum of the 
electron is not a well-defined quantity. However, since U{x) varies smoothly over a 
length scale d ^ Ax, it is possible to introduce a well-defined momentum over the 
size of the small segment under consideration. The expression for the momentum thus 




AF 



(3) 
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depends on the position x of the small segment, and for a state of energy e is given by 
Pe{x) = ±^2m [e — U (x)]. Here the + sign corresponds to the right branch, the — sign 
to the left one. 

Similarly, the eigenstates of the free Hamiltonian are no longer given by simple 
plane waves but instead satisfy —h^d1/2m + U{x) '^e{x.) = e^e{x)- Keeping in mind 
that the typical length scale d associated with the inhomogeneities of the wire is much 
larger than the Fermi wavelength, we use the semiclassical approximation, which yields 

1 



^^\vjx)\ 



exp dx'^2m [e-U{x')] 



(4) 



normalized according to / dx'^ ^^±{x)'^*, 



± 



X) 



27r5(e — e'). Here the velocity is defined 



as v^{x) = p^{x)/m and the index ± refers to the right/left branches. We ignored the 
backscattered wave, since it only leads to exponentially small contributions for kpd ^ 1. 

The rate of change of momentum associated with the three processes shown in 
figure [2] is evaluated using the Fermi golden rule, so that the damping force acting on 
the electrons takes the form 



AF 



2n 

p,k,p',k' 



pk;p'k' 



+ 



y2± 
^ pk;p'k' 



+ 



^ pk;p'k' 



6{ep + efc — ep' — 



x{p' + k'-p-k) [fp^ftil - fp^){l - ft) - f?fUl - - fk)\ , (5) 

where we introduced Vpk^p'k' as the matrix element of the interacting Hamiltonian fl2cj) 
for scattering from the initial state {p, k) to the final state {p', k') according to the 
processes shown in figure [2l The superscripts ||, 2± and 1± refers to the standard 
notations for these scattering processes [35]. The occupation numbers f^'^ introduced 
in ^ are given by the Fermi distribution evaluated with the appropriate temperatures 
Tr^l{x), defined in section [2l 

One readily sees from that the damping force vanishes at Tr = T^. Using the 
fact that the temperature difference Tji — Tl oc I is small in the linear response regime, 
we expand the occupation numbers f^'^ to first order in Tr — Tl. To avoid redundant 
derivations, let us focus on the first process shown on figure EJ^a) . The damping force 
corresponding to this scattering process is then given by 



AFii 



327r^eeF(a;) 



depdekdepidew 



^^pi f fc) ^p' ; ^k' 



2 Epi — ep 



T 



X (p' + A;' - p - A:) 5(e, + - ep. - e,,) - - f^^ 



ipJk\^ Jp'J\^ Jk'Ji (6) 

where we converted the summations over states into energy integrals, and introduced 
the matrix element l^" (cp, e^; ep/, efc/) evaluated using the set of eigenstates defined in 

(H- 

Note that the expansion in — Ti^ leads to an expression for the damping force 
which is proportional to the applied current. As a result, in the linear response regime 
we can ignore any further dependence on J, as this would lead to contributions that are 
non-linear in the current bias. This allows us to use the Fermi energy ep and velocity 
vp they are defined in the equilibrium state, i.e. in the reference frame where the 
electric current vanishes. 
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Focusing on states close to the Fermi energy, we can simplify the expression for the 
eigenstates (jl]) of the free Hamiltonian into 

dx' 



,^,±|^xjexp 



±i(e - e^) / 
Jo 



where '$ep,±{x) is obtained from (jlj) by setting e 



(7) 

hVF{x')\ ' 

ep- This allows us to estimate the 



matrix element to first order in the interaction: 



^" {^p, ^k] V' ^k') 



x+Ax 



dR 



exp 




— + e/c — Cfc' 
hvpix') 



X 



(8) 



[hVF{R)] 

dr V{r, R) (l - 

J -Ax ^ 

Here we introduced the Fermi wave vector kpiR) = pF{R)/k- 

Using the fact that p and p' on the one hand, and, k and k' on the other hand, 
are on the same branch, we express the momentum difference in terms of a difference 
in energy by introducing the density of states. We then define e = ep' — ep + — e^/, 
and perform the remaining energy integrals. Combining the resulting expression for the 
damping force with ([3]), and substituting the matrix element (ED, we obtain the following 
expression for the resistivity associated with the scattering process of figure [2]^a) 



x] 



Q4e'^eF{x)vF{x)n{x)Ax 



-+^-^j^VoiRi)-V2kARi] 



TchVF{Rl) 



X 



0.112- 



TThVF{R2] 



de 



e/AT 



sinh {e/AT) 



exp le 



R2 dx' 



Ri hVpix') 



(9) 



hvp{Ri)hvp{R2) 

The shortened notations Vq and V2kp correspond to the zero-momentum and 2k p Fourier 
components of the potential K(r, R) with respect to its first variable r defined as: 

^i2kp{R)r 



Vo{R) = Jdr V{r,R) and V2kp{R) = / dr V{r,R)e'- 



(10) 



At this stage, it is convenient to rewrite the energy integral in by replacing 
£^ with a second derivative of the exponential term with respect to Ri and i?2, along 
with the appropriate factors of hvp. Performing an integration by parts in the position 
variables leaves us with an expression involving single derivatives of the dimensionless 
parameters Vq{R) / [nhvp{R)] and V2kp{R)/ [T^fi-Vp^R)]. The remaining integral over e 
can be easily simplified by noticing that it is the Fourier transform of a rapidly decaying 
function which only extends over a range of energy comparable to temperature. For 
temperatures T ^ hvp / d, it reduces to a delta function in Ri — R2 which allows us to 
simplify ([9]) to 

1 2 



T 



QAe^ n{x)ep{x) 



Vo{x)-V2kpix) 



(11) 



7lhvp{x) 

where we expanded the remaining position integral to first order in As. 

The contributions corresponding to the remaining two scattering processes can be 
computed following the same steps and are readily obtained from (11 II) by replacing 
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Vo(x) — V2kp{x) with Vo(x) for P2±{x), and with V2fc^(a;) for pi±{x). Combining the 
contributions from all three processes, the final expression for the resistivity in the 
regime of temperatures T ^ hvp/d takes the form 



This expression clearly stresses that the meaningful inhomogeneous quantity is not just 
the interaction potential but rather the dimensionless parameter that involves both the 
electron-electron interaction and the Fermi velocity. In particular, this means that a 
system with a non-uniform density but homogeneous interactions between electrons still 
displays a non-zero resistivity. 

4. Weakly interacting electrons in the moving frame 

We now introduce a different approach for evaluating the resistivity of the system. 
Unlike the derivation of the previous section, this new treatment is compatible with the 
bosonization formalism. Along with providing an alternative derivation of the result 
(fT2|) . our goal in developing this approach is to compare with the results of previous 
works on the inhomogeneous Tomonaga-Luttinger liquid flE l fT9l [20] . 

4.I. Bosonization 

Previous attempts at studying the transport properties of quantum wires relied on an 
extension of the Tomonaga-Luttinger model to inhomogeneous systems fl8\ [T9l [20] . 
These authors assumed that the inhomogeneities do not change the form of the 
Hamiltonian, and can be accounted for by introducing position-dependent velocities 
and Luttinger-liquid parameters. In the general case, however, a rigorous derivation of 
the bosonized Hamiltonian for these systems is still lacking. Here we show how such a 
bosonized Hamiltonian can be derived explicitly in the case of a non-uniform system of 
weakly interacting electrons. 

The standard bosonization formula for weakly interacting fermions involves the 
Fermi momentum as well as a momentum cutoff (see e.g. [IT]), and as such cannot 
be straightforwardly extended to non-uniform systems where both these quantities 
can develop a position dependence. The key idea then is to map the inhomogeneous 
system of electrons onto a set of fictitious fermions described by a Hamiltonian whose 
non- interacting part is translationally invariant. From there, a standard bosonization 
procedure holds and the resulting Hamiltonian expressed in terms of the new variables 
is very reminiscent of the conjectured inhomogeneous Luttinger liquid Hamiltonian, in 
the limit of weak interactions (see section [5]). 

Our starting point is similar to the one we considered in section [3l namely a system 
of interacting electrons with a non-uniform density described by the Hamiltonian (2). 




(12) 



Electronic transport in inhomogeneous quantum wires 



10 



As we noticed in the previous section, the non-uniform potential U{x) appearing in (|2B 
breaks the translational invariance of the system, already when no electron-electron 
interaction is present. As a result, the eigenstates of the free Hamiltonian are no longer 
plane waves but for energies close to the Fermi level, they can be approximated by ([7]). 

Up to a prefactor which depends on position but not on e, these low-energy 
eigenstates look like plane waves, putting forth a more natural set of variables: the 
energy difference e — ep and X{x) = dx' /hvpix'). An expansion of the electron field 
operator ip'y^±{x) over these plane waves calls for the introduction of a fictitious fermion 
field operator rj^^±{X) defined as 

i)^,±{x) = "^epA^) rii,± , (13) 

where '$ep,±{x) was introduced in ([71). Note that the anti-commutation relations satisfied 
by 'i/j^^±{x) transfer to ri^^±{X) ensuring that ^ri^^fj{X),rij^^^,{Y)^ = dcra'^^yisSiX — Y). 

Let us now derive the Hamiltonian describing the physics of these fictitious fermions. 
This is accomplished by substituting (fT3l) into the Hamiltonian (2). By construction, the 
free Hamiltonian is translationally invariant in the new variable X. At low energy, the 
interacting part of the Hamiltonian can be decomposed in three sectors corresponding 
to the conventional gi, g2 and g^ processes [35]. The main difference here is that the 
associated coupling constants are now position-dependent. They can be obtained from 
the Fourier components of the electron-electron interaction potential. 

As an example, consider the so-called g2\\ process. Following [17], the coupling 
constant for this process is given by the zero-momentum Fourier component of the 
interaction potential, which in the case of our inhomogeneous system corresponds to 
Vo(x), introduced in ( ITOl) . Replacing ip with rj according to ( fT3l) . and introducing the 
density operator i>^^j-{X) = ril^^^{X)ri^^±{X), the g2\\ process retains the same form 

dx g2\\ix)p^^^{x)p^_^{x) — >n J dX y2\\iX)uy^^{X)u^_^{X) (14) 

only with a dimensionless coupling constant given by y2\\ (X(x)) = g2\\{x)/7ihvF{x). A 
similar treatment can be applied to the remaining sectors of the interaction. 

The resulting Hamiltonian expressed in terms of the fictitious field rj can now be 
bosonized following the standard procedure: 

V1,±{X) = ^^exp 1^ [±MX) - e,{X) ± MX) - ^.(X)]| (ISa) 

vi,±{x) = ^^exp 1^ [±0,(x) - e,{x) T MX) + e,{x)]^ (155) 

where we introduced the fields (p^, and 61, (with u = p,cr) satisfying bosonic commutation 
relations [(f)^{X), dyOuiY)] = m6{X — Y). Here U^,± are the standard Klein factors [IT] 
and is an energy cutofl[|| introduced to regularize the theory in the ultra-violet sector. 

I Considering that the fietitious fermions were introduced in the vicinity of the electron Fermi surface, 
one should assume a^^ <C ep- 
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In terms of the bosonic variables, the Hamiltonian of the system can be written 
as a sum of two terms describing the excitations of charge and spin degrees of freedom 
respectively, and takes the form 

(16a) 
(166) 



H 


= Hp 




Hp 


1 


fax 




~ 2^ 






1 


/dX 


~ 2^ 





(dxe,f + il + y,iX))idx<l)^f 



+ 



{27iay 



idx9af + (l - vAX)) idx(p.f 
I dXy^iX)cos (2v^0,). 



:i6c) 



The dimensionless parameters yp and y„ are conventional notations for combinations of 
yi, y2 and y^ [IT] given by 

^2fc^(x)-2K,(a;) V2kA^) 



ypixix)) 



(17) 



where Vq and V2kp are the Fourier components of the interaction potential as defined in 

m- 

Note that this form of bosonization clearly highlights that the important variables 
are the dimensionless parameters yp^^ rather than the various interaction constants 5'i,2,4- 
The calculation of the resistivity carried out in section [3] led to the same observation, 
see (H. 



4-2. Resistance and dissipation 

A way to determine the resistivity of the system is to relate it to the mechanism of 
dissipation of energy into the wire when an external current bias is applied. This 
relation was explored in [36] in the context of a quantum wire in the Wigner crystal 
regime, and the method we outline here is similar. 

In the presence of an applied current I = Iq cos ut, the electrons start moving in 
the wire. More specifically, in the dc limit — > 0, one can assume that the current is 
uniform throughout the wire and all electrons move in phase. As a result, the position 
of the electrons depends on time and is related to the injected charge q(t) = IoUJ~^ sin cut 
defined as I{t) = q{t). This time dependence of the positions of the electrons can be 
accounted for by replacing x — x + q(t)/en{x) in the position-dependent parameters of 
the Hamiltonian. While this has no effect in practice when the translational invariance 
holds, for an inhomogeneous system it leads to a time-dependent perturbation to the 
Hamiltonian. Alternatively, this amounts to describing the system in the reference frame 
moving with the electron fluid. In this case, the electrons experience the effect of an 
inhomogeneous potential moving as a function of time. 

In terms of the fictitious set of fermions ri^^±{X), one needs to substitute X in 
the dimensionless interaction parameters ?/p^o- by the time-dependent position X + 
q{t) / en{X) where the density in these variables is given by n{X{x)) = n{x)hvF{x). 
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In the linear response regime, an expansion to first order in q{t) leads to the following 
form of the Hamiltonian: 

en{X) 

where we introduced the notation 7i'(X) = T-C'p{X) +T-C'^{X) for the following quantities: 



J en A 



KiX) = — [dxy.iX)] idx<P,f (19a) 

KiX) = [dxyaiX)] {dx<Paf + TTT^ [dxyaiX)] COs(2v^0.). {19b) 

The time- dependent perturbation in ( |T8|) acts as an external driving force, resulting 
in the creation of spin and charge excitations. These excitations are responsible for 
dissipating the energy from the external force into the wire. Using the Fermi golden 
rule, it is possible to estimate the rates of these absorption and emission processes, and 
therefore, the energy W dissipated in unit time into the system. In the linear response 
regime, where the amplitude Iq of the current oscillations is weak, the energy dissipated 
in unit time is quadratic in Jq and is given by 

h \2euj) J~oo 2-nh ^ ^ 

where (...) stands for thermodynamic averaging. 

The resistance of the system is then derived by comparing the dissipated energy 
obtained in fl20|) with the Joule heat law W = IqR/2. Since the charge part fllQap 
and the spin part (\19b\i of the time-dependent perturbation commute, one expects the 
resistance R of the wire to be expressed as the sum of a spin and a charge contribution 
R = Rp + Ra, which can be evaluated separately. This can be understood as the 
consequence of having two independent channels for dissipating energy throughout the 
wire, corresponding to spin and charge excitations [36] . After some manipulations, these 
two contributions to the resistance can be expressed in the dc limit as 

he^ J n[X) J n[Y j uj 
Here we introduced the retarded correlator yVret,u as the Fourier transform in time of 
WretAX,Y;t) = -ie{t){[K{X,t),nl{Y,0)]), where u = p,a. 



4-3. Charge contribution to the resistivity 

In order to derive the contribution to the resistivity from charge degrees of freedom, we 
substitute in fl2Tl) the expression for Hp introduced in fll9al) . The retarded correlator 
y^ict,p resulting from this substitution is quartic in the bosonic field (pp. It is thus more 
convenient to express it in terms of the corresponding time-ordered correlation function. 
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Wret.p (X, Y; u) = -— [dxy,{X)] [dvypiY)] 



via an analytic continuation in frequency space. This allows for the use of Wick's 
theorem, ultimately leading to the following expression of the retarded correlator: 

n_ 

2" 

(22) 

where {T . . .) corresponds to the time-ordered correlation function. 

Since (1221) is explicitly quadratic in the interaction, the dominant contribution to 
the retarded correlator Wret,p can be derived by using the free propagator of the bosonic 
field (pp. The latter is readily obtained from the non-interacting Hamiltonian, and is 
given by 

(T0,(X,r)0,(F,O)) = TY.J (23) 

Combining (!23|) with (122!) . performing the integral over imaginary time r, and 
substituting the analytically continued result into the charge contribution to the 
resistance in the dc limit reads: 

8e^ J ^(-^) 




1 f 27iT{X - Y) 



(24) 



sinh' [27tT{X - Y)] [ tanh [2nT{X - F)] J ' 

where we restricted ourselves to contributions up to second order in the interaction. 

One recognizes in fl2^ a rapidly decaying integral kernel for |X — y| ^ 1/T. In 
terms of real space quantities, this corresponds to distances of order Tivp/T. It follows 
that at temperatures T ^ hvp/d the double integral in X and Y is dominated by short- 
range contributions. This allows us to reduce the expression for the resistance to 
a single integral over X. Changing variables back from X to x, the integrand of the 
resulting expression for the resistance Rp can be identified with the charge contribution 
to the resistivity at position x in space and is given by 

h T 

''^^^^ = 128e^ eAxMx) ^^^y^^^^^' ^ (25) 
where we focused on temperatures in the range hvp/d ^ T ^ ei?. 



4-4- Spin contribution to the resistivity 

The method used above to derive the charge contribution to the resistivity can be readily 
extended to evaluate that of spin degrees of freedom. Substituting the expression for 
into fl2T]) . one notices that the spin contribution splits off into two parts: one coming 
from the quadratic term in (per, the other from the cosine term. 
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4.4- 1- Quadratic term The expressions for the quadratic parts of H'^ and H'^ are 
identical up to a sign, upon replacing the charge parameter yp and field (pp by their 
spin counterparts. As a result, in order to derive the contribution to the transport 
properties from the quadratic term in 00-, it is sufficient to repeat the steps leading to 
( I25l) but replace i/p by y„ so that 

Pa,quad(a;) = ^ ( \ ( \ [^xVaix)]^ ■ (26) 

128e^ ei?(x)n(x) 

Here we again restricted ourselves to temperatures in the range hvp/d <^ T <^ ep. 

4-4-2- Cosine term The contribution to the resistivity coming from the cosine term 
of the spin Hamiltonian can be inferred from fl26|) based on the following symmetry 
argument. In terms of the bosonized Hamiltonian, the interaction-dependent term 
appearing in the quadratic part of accounts for the coupling between z components 
of the electron spins [IT]. This term ultimately leads to the contribution Po-,quad obtained 
in fl26l) . On the other hand, the cosine term in fll6cp corresponds to the coupling of the 
remaining x and y components [T7]. Because of the SU(2) symmetry, the contributions 
from all three components of the interaction between electron spins are the same. As a 
result, we expect the cosine term to contribute twice as much to the resistivity as the 
quadratic part of the spin Hamiltonian. 

This can be verified explicitly by substituting the cosine term from fllQal) into the 
expression for the retarded correlator Wret.o-- After performing the analytic continuation 
and taking the dc limit — 0, the cosine-cosine correlation function takes the form 



lim — Im 

uj—*0 ijj 



/ dre'^""(Tcos(2V20,(X, r)^ cos (2V20<,(r,O))) 
n^UT'^a^ \ 27rT(X - Y) 



sinh^ \2txT{X -Y)\X tanh [27rT(X -Y)\]' ^^^^ 
Here it is sufficient to perform the thermodynamic averaging using the free Hamiltonian 
since this term enters Wrct.o- with a prefactor quadratic in the interaction parameter y„. 

One recognizes in ( 1271) the same short-range kernel we encountered in (l24l) . It follows 
that, at temperatures T 3> Tivp/d, one can simplify the expression for the resistance 
into a single integral over X. Changing back variables from X to x, and identifying 
the integrand in x with the resistivity, we obtain the contribution from the cosine term 
pa,cos{.x) = 2po-,quad(a^)5 as wc argucd from the SU(2) symmetry. The total contribution 
from spin degrees of freedom thus amounts to three times the result ( 1261) . 

Combining the charge and spin contributions, the final expression for the resistivity 
of the wire at temperatures T ^ hvp/d is given by 

Using flTTl) to replace the dimensionless parameters yp{x) and ya{x) with their expression 
in terms of the electron-electron interaction potential, the latter result becomes identical 
to (da). 
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5. Discussion 

5.1. Inhomogeneous Luttinger liquid 

It is interesting to compare the bosonized Hamiltonian we derived in (16) to that of the 
inhomogeneous Tomonaga-Luttinger model conjectured in [THl [191 EO]. To do so, we 
change variables back from X to x, so that the Hamiltonian (16) takes the form 

H =Hp + H^ (29a) 
\d,i3,f + {l + yp{x)){d,v,f] (296) 



H,= Idx ^"^^^^ 



27C 

H^= fdx \{d,^^f + (1 - y^{x)) (d,^^) 



where we denoted (pu{x) = 0j^(X(a;)) and = 9^{X{x)) (with z/ = p,cr). 

The charge Hamiltonian Hp and the quadratic part of H^j are identical to the 
inhomogeneous Tomonaga-Luttinger Hamiltonian [THl UHl EO], taken in the limit of 
weak interactions. The important difference comes from the cosine term of the spin 
Hamiltonian. This term was absent from previous works which either discarded it 
arguing that the coupling constant g^j renormalizes towards zero at low energy scales 
[37] , or simply focused on a system of spinless fermions [iHl UHl [20] . To recover a standard 
form for the cosine term, one needs to introduce a position-dependent momentum cutoff 
[a(x)]~^ defined as a{x) = ahvpix), where is the energy cutoff introduced in (15). 
Keeping in mind that momentum is no longer a conserved quantity in our model, this 
position dependence of the momentum cutoff was to be expected. Interestingly though, 
the natural guess relying on the common interpretation of a{x) as a small distance 
cutoff, would have led to a different answer. Indeed, assuming that a{x) represents the 
shortest inter-particle distance, one would expect it to be inversely proportional to the 
electron density, i.e. a{x) oc 1/vf{x). 

Although the quadratic part of our bosonized Hamiltonian is similar to the model 
considered [HI [IHl [20], we do not reach the same final answer. This is because our 
treatment amounts to considering perturbations to the Luttinger-liquid Hamiltonian 
which were not taken into account in previous studies. In section [H we treated I{t) 
as an external parameter. In the framework of the Luttinger liquid model, it can also 
be interpreted as an excitation of the charge mode. Using the bosonization expression 
for the electric current I = e{\/2/n)(f)p, one readily sees that q{t) can thus appear as a 
dynamical variable, directly proportional to the charge field (pp. As a result, the linear in 
q{t) perturbation to the Hamiltonian in ( ITSl) corresponds, in the conventional Luttinger- 
liquid theory, to cubic terms in the bosonic fields of the form 4'p{dx4'i>)'^ = p, o"). 
These terms are irrelevant perturbations to the Luttinger-liquid Hamiltonian, and as 
such are usually discarded. However, it was proven that within the quadratic Luttinger- 
liquid Hamiltonian, non-uniform electron-electron interactions do not contribute to the 
resistance [iHl UHl [20] . It thus makes sense to take these irrelevant perturbations into 
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account. Our approach showed that they affect the transport properties in a non-trivial 
way and lead to a finite resistivity. 

5.2. Connection with experiments 

Our results are relevant to experiments performed on long quantum wires. However, 
we focused on the case of weakly interacting electrons which is unlikely to be realized 
in experimental situations. Let us discuss to what extent our conclusions are modified 
when this restriction on the interaction strength is relaxed. 

Though our results are not readily applicable in the case of strong electron-electron 
interactions, the method developed in section H] which relies on bosonization suggests 
that the temperature and density dependences of the resistivity should not be affected 
by the strength of the interactions. Experimental measurements of these dependences 
may thus be compared with our results. 

Furthermore, given the Hamiltonian of the system in the strongly interacting 
regime, one could repeat the treatment of section H] in order to derive the resistivity. 
Unfortunately, a rigorous derivation of the bosonized Hamiltonian in the case of a 
strongly interacting inhomogeneous system is yet to be found. 

5.3. Equilibration 

In our derivation, we assumed that the electronic subsystem is in equilibrium in the 
moving frame. For this to be satisfied, we need the wires to be longer than the typical 
length scale /eq associated with the processes of equilibration taking place inside the wire. 
If the size of the wire becomes too short with respect to the equilibration length /gq, we 
expect our results to be modified by an additional small prefactor of the order of the 
ratio of these two length scales. This might lead to a non-trivial temperature dependence 
of the resistance, depending on the leading equilibration mechanism involved. 

Little is known about equilibration mechanisms in one-dimensional interacting 
systems. In the case of weakly interacting electrons, recent work [M] suggests that 
scattering processes involving three electrons may be the leading source of equilibration 
in the system. Because of consideration of energy and momentum conservation, these 
three-particle collisions should involve states near the bottom of the band, resulting in 
a strong suppression at low temperatures. In the experimentally relevant case of low 
electron density and strong interactions, this analysis no longer holds and a detailed 
treatment remains elusive. It is natural to expect that the equilibration in the wire 
would become easier as the interactions grow stronger. 

6. Summary 

In this paper we studied the effect of inhomogeneous electron-electron interactions on 
the transport properties of a quantum wire. We considered a very general form of 
the interaction potential, and allowed for a non-uniform density of electrons along 
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the wire. We argued that the inhomogeneities allow for non-momentum-conserving 
scattering processes which give rise to a finite resistivity of the wire. We showed that in 
the regime of weakly interacting electrons, such scattering processes contribute to the 
resistivity as a linear in T term|§|, over a broad range of temperatures T below the Fermi 
energy. We also reformulated our results within the framework of the inhomogeneous 
Tomonaga-Luttinger model, and analyzed the differences with previous works relying 
on this formalism. 
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